An improved method for determining frequency of multiple variants of SARS-CoV-2 in wastewater using qPCR assays.

Wastewater-based surveillance has become an effective tool around the globe for indirect monitoring of COVID-19 in communities. Variants of Concern (VOCs) have been detected in wastewater by use of reverse transcription polymerase chain reaction (RT-PCR) or whole genome sequencing (WGS). Rapid, reliable RT-PCR assays continue to be needed to determine the relative frequencies of VOCs and sub-lineages in wastewater-based surveillance programs. The presence of multiple mutations in a single region of the N-gene allowed for the design of a single amplicon, multiple probe assay, that can distinguish among several VOCs in wastewater RNA extracts. This approach which multiplexes probes designed to target mutations associated with specific VOC's along with an intra-amplicon universal probe (non-mutated region) was validated in singleplex and multiplex. The prevalence of each mutation (i.e. VOC) is estimated by comparing the abundance of the targeted mutation with a non-mutated and highly conserved region within the same amplicon. This is advantageous for the accurate and rapid estimation of variant frequencies in wastewater. The N200 assay was applied to monitor frequencies of VOCs in wastewater extracts from several communities in Ontario, Canada in near real time from November 28, 2021 to January 4, 2022. This includes the period of the rapid replacement of the Delta variant with the introduction of the Omicron variant in these Ontario communities in early December 2021. The frequency estimates using this assay were highly reflective of clinical WGS estimates for the same communities. This style of qPCR assay, which simultaneously measures signal from a non-mutated comparator probe and multiple mutation-specific probes contained within a single qPCR amplicon, can be applied to future assay development for rapid and accurate estimations of variant frequencies.

Realizing the value in "non-standard” parts of the qPCR standard curve by integrating fundamentals of quantitative microbiology

The real-time polymerase chain reaction (PCR), commonly known as quantitative PCR (qPCR), is increasingly common in environmental microbiology applications. During the COVID-19 pandemic, qPCR combined with reverse transcription (RT-qPCR) has been used to detect and quantify SARS-CoV-2 in clinical diagnoses and wastewater monitoring of local trends. Estimation of concentrations using qPCR often features a log-linear standard curve model calibrating quantification cycle (Cq) values obtained from underlying fluorescence measurements to standard concentrations. This process works well at high concentrations within a linear dynamic range but has diminishing reliability at low concentrations because it cannot explain "non-standard” data such as Cq values reflecting increasing variability at low concentrations or non-detects that do not yield Cq values at all. Here, fundamental probabilistic modeling concepts from classical quantitative microbiology were integrated into standard curve modeling approaches by reflecting well-understood mechanisms for random error in microbial data. This work showed that data diverging from the log-linear regression model at low concentrations as well as non-detects can be seamlessly integrated into enhanced standard curve analysis. The newly developed model provides improved representation of standard curve data at low concentrations while converging asymptotically upon conventional log-linear regression at high concentrations and adding no fitting parameters. Such modeling facilitates exploration of the effects of various random error mechanisms in experiments generating standard curve data, enables quantification of uncertainty in standard curve parameters, and is an important step toward quantifying uncertainty in qPCR-based concentration estimates. Improving understanding of the random error in qPCR data and standard curve modeling is especially important when low concentrations are of particular interest and inappropriate analysis can unduly affect interpretation, conclusions regarding lab performance, reported concentration estimates, and associated decision-making.

Realizing the value in "non-standard" parts of the qPCR standard curve by integrating fundamentals of quantitative microbiology.

The real-time polymerase chain reaction (PCR), commonly known as quantitative PCR (qPCR), is increasingly common in environmental microbiology applications. During the COVID-19 pandemic, qPCR combined with reverse transcription (RT-qPCR) has been used to detect and quantify SARS-CoV-2 in clinical diagnoses and wastewater monitoring of local trends. Estimation of concentrations using qPCR often features a log-linear standard curve model calibrating quantification cycle (Cq) values obtained from underlying fluorescence measurements to standard concentrations. This process works well at high concentrations within a linear dynamic range but has diminishing reliability at low concentrations because it cannot explain "non-standard" data such as Cq values reflecting increasing variability at low concentrations or non-detects that do not yield Cq values at all. Here, fundamental probabilistic modeling concepts from classical quantitative microbiology were integrated into standard curve modeling approaches by reflecting well-understood mechanisms for random error in microbial data. This work showed that data diverging from the log-linear regression model at low concentrations as well as non-detects can be seamlessly integrated into enhanced standard curve analysis. The newly developed model provides improved representation of standard curve data at low concentrations while converging asymptotically upon conventional log-linear regression at high concentrations and adding no fitting parameters. Such modeling facilitates exploration of the effects of various random error mechanisms in experiments generating standard curve data, enables quantification of uncertainty in standard curve parameters, and is an important step toward quantifying uncertainty in qPCR-based concentration estimates. Improving understanding of the random error in qPCR data and standard curve modeling is especially important when low concentrations are of particular interest and inappropriate analysis can unduly affect interpretation, conclusions regarding lab performance, reported concentration estimates, and associated decision-making.

Statistical framework to support the epidemiological interpretation of SARS-CoV-2 concentration in municipal wastewater.

The ribonucleic acid (RNA) of the severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) is detectable in municipal wastewater as infected individuals can shed the virus in their feces. Viral concentration in wastewater can inform the severity of the COVID-19 pandemic but observations can be noisy and sparse and hence hamper the epidemiological interpretation. Motivated by a Canadian nationwide wastewater surveillance data set, unlike previous studies, we propose a novel Bayesian statistical framework based on the theories of functional data analysis to tackle the challenges embedded in the longitudinal wastewater monitoring data. By employing this framework to analyze the large-scale data set from the nationwide wastewater surveillance program covering 15 sampling sites across Canada, we successfully detect the true trends of viral concentration out of noisy and sparsely observed viral concentrations, and accurately forecast the future trajectory of viral concentrations in wastewater. Along with the excellent performance assessment using simulated data, this study shows that the proposed novel framework is a useful statistical tool and has a significant potential in supporting the epidemiological interpretation of noisy viral concentration measurements from wastewater samples in a real-life setting.

A wastewater-based epidemic model for SARS-CoV-2 with application to three Canadian cities.

The COVID-19 pandemic has stimulated wastewater-based surveillance, allowing public health to track the epidemic by monitoring the concentration of the genetic fingerprints of SARS-CoV-2 shed in wastewater by infected individuals. Wastewater-based surveillance for COVID-19 is still in its infancy. In particular, the quantitative link between clinical cases observed through traditional surveillance and the signals from viral concentrations in wastewater is still developing and hampers interpretation of the data and actionable public-health decisions. We present a modelling framework that includes both SARS-CoV-2 transmission at the population level and the fate of SARS-CoV-2 RNA particles in the sewage system after faecal shedding by infected persons in the population. Using our mechanistic representation of the combined clinical/wastewater system, we perform exploratory simulations to quantify the effect of surveillance effectiveness, public-health interventions and vaccination on the discordance between clinical and wastewater signals. We also apply our model to surveillance data from three Canadian cities to provide wastewater-informed estimates for the actual prevalence, the effective reproduction number and incidence forecasts. We find that wastewater-based surveillance, paired with this model, can complement clinical surveillance by supporting the estimation of key epidemiological metrics and hence better triangulate the state of an epidemic using this alternative data source.

RT-qPCR detection of SARS-CoV-2 mutations S 69-70 del, S N501Y and N D3L associated with variants of concern in Canadian wastewater samples.

SARS-CoV-2 variants of concern (VoC) have been increasingly detected in clinical surveillance in Canada and internationally. These VoC are associated with higher transmissibility rates and in some cases, increased mortality. In this work we present a national wastewater survey of the distribution of three SARS-CoV-2 mutations found in the B.1.1.7 (alpha), B.1.351 (beta), and P.1 (gamma) VoC, namely the S-gene 69-70 deletion, N501Y mutation, and N-gene D3L. RT-qPCR allelic discrimination assays were sufficiently sensitive and specific for detection and relative quantitation of SARS-CoV-2 variants in wastewater to allow for rapid population-level screening and surveillance. We tested 261 samples collected from 5 Canadian cities (Vancouver, Edmonton, Toronto, Montreal, and Halifax) and 6 communities in the Northwest Territories from February 16th to March 28th, 2021. VoC were not detected in the Territorial communities, suggesting the absence of VoC SARS-CoV-2 cases in those communities. Percentage of variant remained low throughout the study period in the majority of the sites tested, however the Toronto sites showed a marked increase from ~25% to ~75% over the study period. The results of this study highlight the utility of population level molecular surveillance of SARS-CoV-2 VoC using wastewater. Wastewater monitoring for VoC can be a powerful tool in informing public health responses, including monitoring trends independent of clinical surveillance and providing early warning to communities.

501Y.V2 and 501Y.V3 variants of SARS-CoV-2 lose binding to bamlanivimab in vitro.

The newly emerging variants of SARS-CoV-2 from South Africa (B.1.351/501Y.V2) and Brazil (P.1/501Y.V3) have led to a higher infection rate and reinfection of COVID-19 patients. We found that the mutations K417N, E484K, and N501Y within the receptor-binding domains (RBDs) of the virus could confer ~2-fold higher binding affinity to the human receptor, angiotensin converting enzyme 2 (ACE2), compared to the wildtype RBD. The mutated version of RBD also completely abolishes the binding of bamlanivimab, a therapeutic antibody, in vitro. Detailed analysis shows that the ~10-fold gain of binding affinity between ACE2 and Y501-RBD, which also exits in the high contagious variant B.1.1.7/501Y.V1 from the United Kingdom, is compromised by additional introduction of the K417/N/T mutation. Mutation of E484K leads to the loss of bamlanivimab binding to RBD, although this mutation does not affect the binding between RBD and ACE2.